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Abstract 

HE  velocity  and  vorticity  fields  around  an  oscillating 
airfoil  in  compressible  dynamic  stall  are  reported.  Phase 
averaged,  two  component  laser  velocimetry  data  were  ob¬ 
tained  at  a  freestream  Mach  number  of  0.3  and  a  reduced 
frequency  of  0  05.  This  is  the  first  set  of  velocity  data  available 
at  a  high  Reynolds  number  (540,000)  under  compressible  flow 
conditions  and  it  serves  as  a  good  database  for  development 
and  validation  of  computer  codes.  Of  particular  interest  is  the 
formation  of  a  separation  bubble,  which  bursts  coincidentally 
with  the  formation  of  the  dynamic  stall  vortex,  adding  an 
extra  degree  of  physical  complexity  to  the  problem 

Contents 

The  performance  of  helicopters  and  the  maneuverability  of 
fixed  wing  aircraft  is  limited  by  dynamic  stall  and  the  gross 
separation  of  the  unsteady  flow  that  accompanies  it.  The 
dynamic  stall  flowfield  is  a  complicated  combination  of  a 
multitude  of  fluid  dynamic  effects  that  include  tremendous 
fluid  acceleration  around  the  leading  edg-  formation  of 
strong  suction  peaks;  development  of  the  local  boundary  layer 
under  strongly  adverse  pressure  gradients;  transition  of  the 
leading-edge  laminar  boundary  layer,  its  separation,  and  reat¬ 
tachment  resulting  in  a  separation  bubble  that  grows  and 
eventually  bursts  just  when  the  dynamic  stall  vortex  forms; 
under  compressibility  conditions  (A/  2 0.3),  formation  of 
shock(s)  and  the  induced  separation  due  to  it;  addition  of  large 
amounts  of  coherent  vorticity  into  the  flow  and  its  coalescence 
into  the  dynamic  stall  vortex  that  later  convects  over  the  airfoil 
upper  surface  and  interacts  with  the  trailing  edge  separated 
flow  that  propagates  toward  the  leading  edge;  and  so  on.  The 
complex  interactions  governing  the  physical  problem  have 
hitherto  made  understanding  the  origin  of  dynamic  stall  and 
attempts  at  modeling  the  flow  very  challenging.  The  study 
being  reported  provides  a  comprehensive  set  of  velocity  data 
along  with  vorticity  distributions  derived  from  it,  which  can  be 
used  to  model  the  flow.  In  addition,  it  also  serves  as  a  data¬ 
base  to  verify  computational  results  and  enables  development 
of  new  codes. 


Presented  as  Paper  91-1799  at  the  AIAA  22nd  Fluid  Dynamics, 
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1991;  received  Sept.  20,  1991;  synoptic  received  Sept.  25.  1992;  ac¬ 
cepted  for  publication  Nov.  2,  1992.  Full  paper  available  from  A1AA 
Library,  555  West  57th  St..  New  York,  NY  10019.  Copyright  ©  1991 
by  the  American  Institute  of  Aeronautics  and  Astronautics,  Inc.  No 
copyright  is  asserted  in  the  United  States  under  Title  17,  U.S  Code. 
The  U.S.  Government  has  a  royalty-free  license  to  exercise  all  rights 
under  the  copyright  claimed  herein  for  Governmental  purposes.  All 
other  rights  are  reserved  by  the  copyright  owner. 

•Associate  Director  and  Research  Associate  Professor,  Navy- 
NASA  Joint  Institute  of  Aeronautics,  Department  of  Aeronautics  and 
Astronautics,  Mailing  Address:  M.S.  260-1.  NASA  Ames  Research 
Center,  Moffett  Field.  CA  94035-1000.  Associate  Fellow  A1AA. 

tResearch  Scientist;  on  leave  from  The  National  Aeronautical  Lab¬ 
oratory,  Bangalore,  India. 


A  7.62-cm-chord  NACA  0012  airfoil  was  oscillated  sinu 
soidally  with  its  angle  of  attack  varying  as,  a  =  10  deg- 10  deg 
sin  uil  at  a  reduced  frequency  [A:  =  (ir/c  C'»)l  of  0.05  in  a  flow 
with  a  freestream  Mach  number  of  0.3.  The  details  of  the  wind 
tunnel  used  are  given  in  Ref.  1  A  standard  two  component  ( L 
and  V)  frequency  shifted  laser  velocimeter  was  used  for  the 
measurements.  Optical  encoders  mounted  on  the  drive  pro¬ 
vided  the  pertinent  time-dependent  information.  A  novel  ap¬ 
proach  of  freezing  the  encoder  data  when  coincident  L  and  I' 
velocity  samples  were  present  served  to  provide  reliable  phase 
angle  data,  which  were  later  ensemble  averaged  by  sorting  in'o 
bins  covering  3  deg  in  phase  angle.  The  velocities  were  com¬ 
puted  for  bins  containing  at  least  50  samples.  When  this 
condition  was  not  satisfied,  suitable  interpolation  methods 
were  used,  if  data  were  found  in  the  neighboring  bins.  One  pm 
polystyrene  latex  particles  suspended  in  alcohol  were  used  for 
seeding  the  flow  .  Further  documentation  of  the  experimental 
details  can  be  found  in  Ref.  2. 

Velocity  Measurements  in  the  Separation  Bubble 

Some  interesting  flow  features  can  be  seen  in  the  variation 
of  the  horizontal  component  of  velocity  V  with  phase  angle  o. 
in  Fig.  1,  for  different  heights  (>/c)  at  the  streamwise  location 
x/c  =0.083.  At  y/c  =  0.067,  phase  angle  of  160  deg.  a  =  6.6 
deg,  the  velocity  drops  rapidly  as  a  separation  bubble  forms 
over  the  airfoil.  (It  should  be  noted  that  the  flat  portions  of 
the  distributions  for  y/c  =0.067  and  y/c  =0.083  are  caused 

A  Y/C  =  0.200  O  Y’C  =0.150  K  Y'C  =  0.100 

+  Y/C  =  0.183  VYC  =  0.133  ♦  Y'C  =  0.083 

X  Y'C  =  0.167  B  Y'C  =  0.117  ®  Y'C  =  0.067 
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Fig.l  1/  component  velocity  measurements  at  x/c  =  0.083  (velocity 
offset  by  0.2  (l//l/»>  for  eich  successive  y/c  location). 
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by  the  blockage  of  the  beams  close  to  the  airfoil.)  Since  the 
bubble  is  encountered  later  in  the  cycle  at  a  higher  y/c  point, 
the  phase  angle  at  which  the  drop  occurs  increases  with  y/c. 
However,  at  y/c  values  outside  the  bubble,  the  velocity  distri¬ 
butions  are  parallel  to  each  other.  This  leads  to  the  conclusion 
that  the  bubble  height  is  0.03  c-0.04  c  above  the  airfoil  upper 
surface. 

An  analysis  of  the  corresponding  vertical  component  of 
velocity  V  showed  rapid  increases  ,n  the  velocity  at  <t>  =  200 
deg,  a  =  13.4  deg.  This  is  associated  with  bursting  of  the 
separation  bubble.  It  is  worthwhile  mentioning  that  the  bubble 
bursting  is  somewhat  gradual  and  not  as  abrupt  as  is  normally 
perceived. 
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Fig.  2  '"wtonri  of  normalized  z  component  of  vortidty:  a)  0  —  171 
4eg,  a  -  *.44  deg,  bi  *  =  ISO  deg,  a  =  10.0  deg,  and  c)*  =  198  deg, 
or -  13.09  deg. 


Vorticity  Distributions 

The  2  component  of  vorticity  was  calculated  from  the  mea¬ 
sured  U  and  V  components  of  velocity  by  first  fitting  a  cubic 
spline  curve  to  the  data  and  interpolating  the  velocities  in  a 
grid  at  a  resolution  of  1.25  mm,  using  a  second-order  central 
differencing  scheme.  Thus,  the  noise  level  in  the  distributions 
is  expected  to  be  high  at  about  20%  of  the  local  maximum 
vorticity  values  (in  both  the  positive  and  negative  quantities). 
The  following  discussion  about  the  vorticity  field  should  still 
be  valid,  especially  before  the  dynamic  stall  vortex  begins  to 
convect  because  no  discontinuities  such  as  shocks  were  en¬ 
countered  within  the  measurement  grid.  Thus,  the  picture  of 
the  flowfield  is  also  quantitatively  valid  up  to  the  point  where 
the  particles  were  able  to  follow  the  flow  properly. 

Figure  2a  shows  that  at  d>  =  171  deg  (or  =  8.44  deg),  a  region 
of  clockwise  vorticity  has  developed  over  the  airfoil,  just 
around  the  location  of  the  separation  bubble,  with  a  peak 
vorticity  of  -  8  units  in  it.  A  region  of  counter-clockwise 
vorticity  could  also  be  found  above  it,  but  the  peak  vorticity 
within  it  is  only  about  5  units.  As  the  airfoil  reaches  an  angle 
of  attack  of  10  deg  (Fig.  2b)  the  clockwise  vorticity  has  in¬ 
creased  to  -  1 1  units,  whereas  the  anticlockwise  vorticity  is 
still  at  5  units.  The  extent  of  the  vortical  region  has  grown  to 
about  25%  chord  in  both  the  x  and  y  directions.  As  the  airfoil 
pitches  to  higher  angles  of  attack,  the  vorticity  should  increase 
steadily  until  the  vortex  begins  to  convect.  Figure  2c  shows 
that  at  <t>  =  198  deg,  this  is  the  case  as  the  clockwise  vorticity 
has  doubled  to  -  22  units,  but  the  anticlockwise  vorticity  has 
only  increased  to  about  10  units.  Earlier  experiments’  have 
shown  that  the  vortex  begins  to  convect  at  around  this  phase 
angle.  The  separation  bubble  also  bursts  around  the  same 
angle  of  attack.  Thus,  a  combined  effect  is  felt  by  the  airfoil, 
which  should  be  seen  in  its  vorticity  field.  A  calculation  of  the 
circulation2  over  the  measurement  region  showed  an  increase 
until  the  stall  vortex  convection  was  initiated  and  dropped 
slightly  beyond  this  angle  of  attack. 

The  study  leads  to  the  following  conclusions. 

1 )  One  of  the  salient  features  of  the  flow  is  the  formation  of 
a  separation  bubble.  This  bubble  bursts  (opens  up)  just 
around  the  angle  of  attack  at  which  and  the  location  where  the 
dynamic  stall  vortex  forms,  complicating  the  flow  physics. 

2)  The  clockwise  vorticity  was  found  to  increase  in  the  flow 
until  the  vortex  begins  to  convect. 
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